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ABSTRACT 

We study gaseous outflows from disk galaxies driven by the combined effects of ram pressure on cold 
gas clouds and radiation pressure on dust grains. Taking into account the gravity due to disk, bulge 
and dark matter halo, and assuming continuous star formation in the disk, we show that radiation or 
ram pressure alone is not sufficient to drive escaping winds from disk galaxies, and that both processes 
contribute. We show that in the parameter space of star formation rate (SFR) and rotation speed of 
galaxies, the wind speed in galaxies with rotation speed Vc < 200 km s^^ and SFR < 100 M© yr^^, 
has a larger contribution from ram pressure, and that in high mass galaxies with large SFR, radiation 
from the disk has a greater role in driving galactic winds. The ratio of wind speed to circular speed 



can be approximated as 
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We show that this conclusion is 



borne out by observations of galactic winds at low and high redshift and also of circumgalactic gas. 
We also estimate the mass loading factors under the combined effect of ram and radiation pressure, 
and show that the ratio of mass loss rate to SFR scales roughly as v~^'S~'^, where Eg is the gas column 
density in the disk. 

Subject headings: galaxies: starburst — galaxies: evolution — intergalactic medium 



1. INTRODUCTION 

Galactic winds have been observed at different wave- 
lengths in galaxies of various masses and in a range of 
redshifts. Galaxies, especially with star formation rates 
(Ssi^fl) > 10^^ Mq yr~^ kpc~^, often show large out- 
flow of hot gas that emits X-rays and in which cold clouds 
are found to b e embedded, which a r e observed w ith Ha 
or NaD lines (jHeckman et all 120001 : iMartinI 120051 ). The 
speed of the clouds in the wind range from a few tens to 
several hundred km s~^, and the total mass loss r ate can 
be se veral times the star formation rate ( Veilleu x et al.l 
[20051). 



These outflows play a crucial role in the evolution 
of galaxies by expunging gas, and thereby suppressing 
the star formation. The attempts to understand galac- 
tic evolution in the cosmological context have since long 
encountered the so-called 'cooling catastrophe' problem, 
since left to its own device the baryonic gas would cool 
and form stars more rapidly than observed. It is gen- 
erally believed that a feedback loop inhibits this, and 
that the process of star formation excites an outflow and 
quenches itself. The observed mass-metallicity relation in 
galaxies also indicate that galactic outflows play a major 
role in the chemical evolution in galaxies. Furthermore, 
these outflows enrich the intergalactic medium with met- 
als. 

The standard model to understand galactic outflows 
involves a heated interstellar medium (ISM) under the in- 
fluence of supernovae (SN), and the h ot gas be ing driven 
by th ermal pressure ([Chevalier fc Gleg g 198?^ lHeckman[ 
I2002D . The expansion speed of this hot gas can be large 
enough to eject it out o f the galaxy (|Larsonlll974l : ISaitol 
[19791 iDekelfc Silk|[T986l) . The observ ations of cold gas in 
these outflows ( Heckman et al.|[2000D led to the proposal 
that the cold gas entrained in the hot gas moved due to 
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ram pressure. The wind sp eed was however not found to 
corre late with galaxy mass ([Heckman et al.ll2000t [Martini 
I1999D . and it was argued that the supernovae rate in- 
creased with SFR and hence the wind velocity might 
correl ate with SFR. Simulations also supported this sce - 
nario (iSuchkov et al.l 119941: [Strickland fc Stevens! [2000). 
However, there is a limiting cloud speed implicit in this 
process since ram pressure acts on the cold gas until the 
cold gas velocity becomes equal to that of hot gas. 

This scenario, however, has met with problems from 
new observations of cold component which show that 
the terminal outflow speed depends on g alactic prop- 
erties like rotation speed ([Martini 120051 : iRupke et al.l 
120051) . It has been proposed that these observations can 
be explained by r adiat io n pressure dr i ving the out- 
flow iMurrav et al.l (|2005[ ) : [Marth] (pOOl : [Sha7ma et al.l 
([20111 ). It has also been pointed out that a natural course 
of events leading from a starburst would be a radiation 
pressure driven wind in the beginning, and ram pressure 
acting on it aft er a period of 3-5 My r, the life time o f 
massive stars ([Nath fc Silk 1 120091 : IMurrav et 311120111 ). 
This scenario also naturally explains the puzzling fact 
that cold clouds are observed at large distances although 
their survival time-scales in the hot gas would have in- 
hibited them from being pushed out to such distances. 

In the face of two processes leading to outflows, one 
wonders if both processes contribute equally, or if there 
are regimes in which one of these two processes dominate 
over the other. In this paper we present an analytical 
calculation for the dynamics of cold clouds taking into 
account both ram and radiation pressure and all sources 
of gravity, and compare our results with observations. 

2. GASEOUS OUTFLOWS WITH RAM AND RADIATION 
PRESSURE 

We consider the dynamics of cold clouds (T < 10^ K) 
embedded in hot gas, in which the hot gas component 
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exerts a drag force due to ram pressure. We also assume 
that dust grains in the cold clouds are strongly coupled to 
the gas, and therefore the dynamics of these clouds is also 
influenced by radiation pressure. We ignore magnetic 
forces and the compression of cold clouds by hot wind 
gas. 




Fig. 1. — A schematic diagram for the motion of a cold cloud cm- 
bedded in a cone of hot flow and acted upon by forces of radiation 
and gravity from the parent galaxy. Cloud is at a height z. The 
total mass loss in hot flow, A/^ = 2 p^VftA, where A is the area at 
the top of conical patch, and a factor of 2 for two-sided mass loss. 
For a half-cone angle of 9 the area A = f^^ J^sin{9')d9'd4i. 

We therefore solve the following equation for the mo- 
mentum of the cold cloud (Pc = Mcv), see figure [U 



dPr dv 

£ — M — p 

^^-t^c J, ^ ram 

at at 



grav 



(1) 



where Mc is the mass of the cloud and v is its velocity in 
z direction. Fram represents the force exerted by the hot 
wind via ram pressure in g cm s~^. Frad is the force due 
to radiation on dust grains and Fgrav is the gravitational 
force. 

We first discuss the role of ram pressure on the motion 
of cold blobs of gas dragged in it, following the model of 
iStrel'nitskii fc SunvaevI (|1973[ ). In this scenario, the hot 
gas observed in X-rays and which is thought to provide 
the r am pressure, ha.s temperatur es in the range 0.5-1 
keV ()Heckmanll2002l : lMartinll2005[ ). which correspond to 
the isothermal sound speed Cg ~ 300-400 km/s. Current 
X-ray instruments can not detect the speed of this hot 
and tenuous material and hence the kinemetics of this 
hot phase is poorly constrained. If we assume it as an 
adiabatic wind passing through a sonic point, then vf^ ^ 
acl, where a = 2.5-5 (jEfstathioull2000l ) . which gives Vh ~ 
1.2-2.2 Cg- In this paper, we take Vh ^ 800 km s"-'^, which 
corresponds to Vh ~ 2cs and Tx ~ 1 keV. 



Consider then the hot gas flow (with density ph and 
velocity Vh), emerging through a cone. Mass loss in a 
hot wind is given by the following expression (sec figure 



Mh = 2phVhZ^ / / sin{e')de' dc^ (2) 
Jo Jo 

Observations indicate conical 

in the range 26 ^ 10 °-100° , 

iLehnert fc Heckmanlll996f ). We consider a mass loss rate 
of « {n/2)z^PhVh, which roughly corresponds to half- 
cone angle ^ 30°. The momentum injection rate is 
■ph = MhVh, so we can write 



a ngles for hot wind 
(|Veilleux et all [20051: 



Phvl 



Ph 



7rz2/2 



(3) 



The force exerted by the ram pressure on a cold cloud of 
mass Mc and cross-section Ac is given by 



Here Tiivh 
1 for V < 



^CoAcPhivh 



1 — nivh 



v) (4) 



- v) is the step function whose value is 
Vfi and otherwise. Co ^ 0.5 is the 



drag coefficient. For the cloud, one can write = 
p,mpNH where N// is the column density and p is the 
mean molecular weight. Also the momentum injection 



rate ph is 



5 X 10 
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(t 



SFR 



(jLeitherer et all Il999t) 
eqn.Q in eqn.Q we get. 



Mg/yr J\ dyne in a starburst 
Using these and substituting 



Mc 



[5xlO=^^(T^)]dyne 
4 Nh unip (7rz2/2) 



1-L{vh - v) 
(5) 

Next we consider the forces due to a galactic disk. We 
will use / for force per unit mass (/ = ■^)- In cylindrical 
geometry, the force of gravitation fg^d{z), and that due 
to radiation fr^d{z), along the pole of a disk of radius r^, 
with constant surface density (S) and surface brightness 
(/) are given by. 



fg,d — 2ttGT. 



zrdr 



(r2 



2)3/; 



J, 



2-1^1 z'^rdr 



r,d 



(r2 + z2)2 

where k is the average opacity of a dust & gas mixture. 
The ratio of these forces, the Eddington ratio, increases 
with the height z, beginning with a value of To = 2cgt, 
at the disk centre at z = 0. Since oc tt-, where Mw is 
the disk mass, we can express Fg in terms of the SFR by 
calculating the luminosity L of a galaxy in any desired 
band for a certain SFR using the Starburst99 code. The 
luminosity in this case is proportional to SFR, therefore 
if Li is the luminosity at 1 Gyr for an SFR of 1 Mq yr~^ 
then we can write Fq as, 



Fn = 



Li X 



SFR 

taTqT^ 



(7) 



2cG Md 

We use the mean opacity for gas mixed with dust ~ 200 
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cm^ g ^ corresponding to a color t empe rature ~ 9000 K 
in the U band (Figure Ib. lDraind (|2011[ ). 

To determine the gravitational force, we assume a 
spherical mass distribution in the bulge and halo. For 
the bulge, we assume a total mass of Mi, ^ O.IAI^ in- 
side a radius ri, O-lr^ for simplicity. For the halo, we 
consider a Nava rro- Frenk- White (NF W) profile, with to- 
tal mass Myir ([Navarro et al.|[T997[ ). We fix Myir for a 
given disk mass (M^), by the ratio Myir/Md ~ 20, as 
determined by Mo, Mao & White (1998) (referred to as 
MMW98 hereafter). We evaluate the disk exponential 
scale-length (r^) using the prescription of MMW98, and 
use it as the size of galactic disk. Gravitational potential 
of NFW halo is. 



GMy 



NFW 



ln(l + c)-c/(l+c) 



ln(l + f) 
R 



(8) 



where R = y/r^ + z^, c ~ 



■^^20 ig the concentration pa- 
rameter, is the NFW scale length and R200 is the 
radius within which the mean overdensity is 200. This 
potential implies a gravitational force along z which is 
given by. 



fhalo,z 



NFW 



dz 



(9) 



ln(l + -^)-z/(z + i?,) 



ln(l + c)- c/(l + c) 



The rotation speed implied by the NFW profile peaks at 
a radius R ^ 2Rs, given by. 



ln(3) - 2/3 



2 ln{l + c)-c/{l + c) 



(10) 



where V200 is the rotation speed at R2oo- We choose this 
value of the maximum rotation speed to represent the Vc 
of the disk galaxy, since Figure 2 of MMW98 shows that 
the value of Vc from the fiat part of the total rotation 
curve does not differ much from the peak of the rotation 
curve from halo only. The escape speed in a NFW halo 
is given by 



ln(3) 



R2 



R 



■ln(l 



cR 



R 



200 



) 



1 + c 



(11) 

Figure [2] shows the escape speed along the z-axis for dif- 
ferent galaxies. The dashed and solid lines show the es- 
cape speed at 10 and 20 kpc from the disk plane, for 
galaxies with different circular speed. We have used the 
relation between the halo concentration param eter c and 
galactic mass as given bv lMaccio et"all (I200I . We find 
that for low mass galaxies with Vc < 100, the escape 
velocity Vf.sc ^ 2uc, and that for higher mass galaxies, 
the escape speed ranges between 2-3 Wc- We can there- 
fore conclude that for escaping winds, the ratio of wind 
speed to circular speed should be in the range of 2-3. 

One can finally rewrite eqn(IT]) for evaluating the veloc- 
ity of clouds as a function of z, 
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Fig. 2. — Ratio of NFW escape speed to the circular speed vs. the 
circular speed, at two different vertical distances, 10 kpc (dashed 
line), and 20 kpc (solid line). Dotted line is for a fixed value of 
halo concentration parameter c = 10. 
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where Fq is given by eg. [7 1 . 

10^1 cm-2 ([Martini l2005rfHeckman et al.l[2000[ ). Here 



ln(l + c)-^ 
We use iJ. = lA and 



the first term on RHS denote ram pressure, second the 
radiation pressure and the last three terms represent the 
gravity of the disk, bulge & NFW halo respectively. This 
equation is non-linear due to the presence of v in ram 
pressure term and should be solved numerically, although 
previous authors have approximated it assuming w <C Vh- 
The form of the ram pressure term suggests that ram 
pressure would not be effective once the velocity becomes 
greater than velocity of hot component. Hence the ram 
pressure is likely to be effective for low-mass galaxies. 
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Fig. 3. — Variation of wind speed with vertical distance {z) for 
galaxies of different circular speeds. The thick solid, dashed and 
dotted lines refer to Vc = 100 km s-\ and for SFR of 10, 100, 500 
Mq yr~^, respectively. The thin solid, dashed and dot-dashed lines 
refer to Vc = 200 km s~^, for the same values of SFR, respectively. 
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3. RESULTS 

We solve the wind equation fean ll2p numerically. Fig- 
ure[3]shows the wind speed as a function of z for different 
values of SFR for two galaxies, with Vc = 100 km s~^ and 
Vf, — 200 km s^^. Instead of rising continuously, the wind 
speed saturates after travelling a distance of > 10 kpc, 
with a terminal speed that is lower for higher mass galax- 
ies. The thick solid line roughly corresponds to M82, and 
the wi nd speed ~ 300 km s"-'^ is consistent with obser va- 
tions (jHeckman et al.ll2000t [Schwartz fc Martinll2004[ ). 
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Fig. 4. — Ratio of wind velocity at 20 kiloparscc and the galactic 
rotation speed plotted with Vc for three different SFR in the left 
panel & with SFR for three different Vc in the right panel. Three 
representative cases arc shown with solid circle (dwarf starbursts), 
solid triangle (LIGs) and solid square (ULlGs). The thin horizontal 
line corresponds to ti^, = 3vc. 

We then use the wind speed at z = 20 kpc and show 
the variation of Vw/vc with circular speed Vc and SFR 
in the left and right panels of Figure H] respectively. We 
find that, for a constant SFR, v^jvc decreases with Vc, 
as gravity increases with Vc- We also show three rep- 
resentative cases in this plot, of dwarf starbursts (solid 
circle: - 60 km s-\ SFR - 5 M© yr'^), LIGs (solid 
triangle: Vc ^ 120 km s"!, SFR ~ 50 Mq yT~^), ULIGs 
(solid square: ~ 220 km s-\ SFR - 500 Mq yr"!). 
The values of v^/vc lie close to 3 which is shown by 
the thin horizontal line. The near constancy of Vw/vc 
for the three representative points recovers the observed 
scaling of v.^, with Vc. Taking into account the variation 
of Vw/'Vc with Vc and SFR we find that, the results can 
be approximated by the following fit. 
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(13) 



Next we solve the wind equation for a grid of SFR and 
galaxy circular speed values, for the cases of ram pressure 
and radiation pressure alone, and then for the combina- 
tion of the two. In Figure [5] the wind velocity is zero in 
top left corner for high mass & low SFR galaxies. Wind 
velocity increases as one moves diagonally, from top left 
to the bottom right corner. We show two contours for 
Vu, = 3vc with thin solid lines for ram and radiation pres- 
sure alone. For the case of combined ram and radiation 
pressure driving, we show two contours, for = Vc and 
ivc (upper and lower thick lines). We also show the data 
for outflows from a number of observations (see caption 
for details). 
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Fig. 5. — Contours of f^_20fcpc in units of Vc for winds driven by 
only ram pressure, only radiation pressure & the combination of 
two. The wind velocity increases as one moves from top left corner 
to bottom right corner. Also plotte d are the data po ints with 
different symbols: hollow c ircles (He ckman et al.l 120001). square s 
IIMartinI 12003; 'Genzcl ct al. 200 f). triangles llWeiner et al.|[2009l '). 
diamonds (Schwartz & Martin 20(3), big f^d cr oss with circle at 
its ce nter (iRupke ct al.n2005i )). big blue cross IITumlinson et ahl 
II2011I ) without two outliers in SFR). 



In the case of only radiation pressure, the wind speed 
is found to be roughly proportional to SFR, which can 
be understood from the fact that Fq oc SFR. The case 
for only ram pressure appears to explain the wind in low 
mass galaxies. However, from the Vw = Stic contour it is 
clear that ram pressure can not drive the cold clouds out 
of the galaxies with rotation speeds > 200 km/s, as we 
have seen in the previous section that for escaping winds 
one needs ^ 3uc. This points to the existence of a 
critical rotation speed above which the physical mecha- 
nism of outflow changes. Therefore outflows from galax- 
ies with Vc < 200 km/s and SFR < 100 Mg/yr are dom- 
inated by ram pressure and those from the more massive 
galaxies with larger SFR, are influenced more by radia- 
tion pressure. 

4. DISCUSSIONS 

The most important result of our calculation is that 
galactic outflows require both ram and radiation pres- 
sure, especially for high mass and high SFR cases. Our 
calculation has a number of ingredients from stellar 
physics and disk and halo parameters, and apart from 
the value of the hot wind speed Vh , there is no free pa- 
rameter in this calculation. It is therefore interesting 
to note that our theoretical results are consistent with 
most data of outflows when studied in the parameter 
space of Vc and SFR. It is also interesting that a re- 
cent simulation with ram and radiation pressure driven 
outflows has concluded that these two processes are im- 
portant in different mass regimes, although it is not 
clear wh ere the dividing line between the two regime s 
lies (iHopkins et al.ll2011t Fron de Voort fc Schavdl2011[ ). 
Cold cloud outflows from galaxies on the left of the con- 
tours in figure [H] are unlikely to escape into the IGM 
and likely get trapped in the circu mgalactic region as ob- 
served by iTumlinso n et al.* (120111) fda ta shown by blue 
cross) or fall back (Oppenhcimer fc Davg, 20081) . 
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Although strictly speaking our calculation refers to 
cold clouds being driven out along the pole of the disk 
galaxies, and we cannot infer the mass loss rate with- 
out doing a 2-D calculation, but we can speculate on 
the scaling of the mass loss rate with galactic mass by 
making some simple assumption. Let us assume that 
the dynamics of cold clouds beyond the polar regions 
are similar to that along the pole. Assuming a one- 
dimensional mass flow, the mass loss rate from the disk 
is approximately oc t;u,[Sg7rr^], where is the gas 
column density and is the scale length of the disk. 
We note that in the prescription of MMW98, one has 
Vc oc rd- We therefore have, Mw oc w^~°'^^M°-^I]g, where 
we have used eqn[l31 after multiplying both sides by Vc- 
The ratio of mass outflow rate to the SFR is therefore 
Mw/M, oc vl-'^'^T,gM^'^-^. Using Kennicutt's law of star 
formation, which gives oc S^-^r^ oc E^'^u^, we have 

finally, oc vl'^^'Eg^-^^ . We can therefore conclude 
that roughly. 



Interestingly, similar power law dependenc e has also been 
found in simulations (jHopkins et al.|[2011|) . 

We note that our results assumed a value of Vh ^ 800 
km s^^, and a column density of cold clouds of ~ 10^^ 
cm^^. If we assume a larger value of Vh (~ 1000 km 



s~^), then the contour for only ram pressure will be able 
to explain the winds in ULIGs with large SFR and high 
mass. A similar result will follow from larger values of k 
for the radiation pressure case. 

It is interesting to note that the contour for only ra- 
diation pressure can explain the ULIG region of Figure 
[5] (top right corner). Extending to larger SFR, our re- 
sults indicate that radiation pressure will also be im- 
portant for HLIGs (Hy perluminous Infra-red galaxies) 
(| Rowan- Robinsonll2000( ) . Lastly, although it may appear 
that the role of radiation pressure in galaxies other than 
ULIGs is less dominant than ram pressure as far as en- 
ergetics is concerned, radiation pressure may still play 
an important role in lifting the clouds to a large height 
befor e it is embedded in the hot wind to hel p it survive 
long (|Nath fc Silk \\2(M . iM^av et al.l[201ll) . 

5. SUMMARY 

We have studied the outflows from disk galaxies driven 
by ram and radiation pressure and compared the theo- 
retical results with data in the parameter space of galaxy 
circular speed and SFR. We found that the driving mech- 
anism of escaping wind is different in low mass and high 
mass galaxies, with radiation pressure being important 
for high mass galaxies with high SFR. Our results are also 
consistent with recently observed circumgalactic gas. 

We thank Mitchell Begelman, Bruce Draine, Tim 
Heckman, Yuri Shchekinov and an anonymous referee 
for valuable comments. 
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